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Simulations from a new neutrino event generator 
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^'Institute of Theoretical Physics, Wroclaw University, 
pi. M. Borna 9, 50-204 Wroclaw, Poland 

We construct a new Monte Carlo generator of events for neutrino interactions. The dynamical models for quasi¬ 
elastic reactions, A excitation and more inelastic events described by the DIS formalism with the PDFs modihed 
according to recent JLab data are used. We describe in detail single pion production channels, which combine 
the A excitation and DIS contribution. Many comparisons of the outcome of simulations with experimental data 
are presented. 


1. INTRODUCTION 

We present a new Monte Carlo generator 
of events for neutrino interactions. The orig¬ 
inal motivation for our work was to improve 
NUX-I-FLUKA scheme where no separate reso¬ 
nance contribution is present The aim of 

NUX-I-FLUKA was to describe interactions of 
neutrinos of higher energies and from that point 
of view the resonance part was of minor impor¬ 
tance. Usually MC generators contain a reso¬ 
nance contribution described by means of Rein- 
Sehgal model covering the kinematical region of 
hadronic invariant mass M + < TU < 2 GeV. 

If for neutrino reactions the quark-hadron dual¬ 
ity holds true one can assume that contributions 
from higher resonances are averaged by deep in¬ 
elastic scattering (DIS) structure functions and 
that only the dominant A resonance has to be 
treated separately. 

The current version of the generator includes 
various dynamical models: quasi-elastic [ 2 ], A ex¬ 
citation P], and DIS for which we use GRV94 Par- 
ton Distribution Functions (PDF) ^ with mod¬ 
ifications proposed by Bodek and Yang [ 5 ]. The 
total cross section for the neutrino scattering is 
assumed to be the incoherent sum: 



where aspp is the sum of cross sections for single 
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pion production (SPP) and CC and NC denote 
charge and neutral current reactions respectively. 

The MC generator is organized around the 
event structure which contains three vectors of 
particles: incoming, temporary and outgoing. It 
also contains a structure with all the parameters 
used and a set of boolean flags tagging the event 
as QEL, DIS, CC, NC etc. The input parame¬ 
ters are read at startup from a text file and the 
events are stored in the ROOT tree file to simplify 
further analysis. Different interactions are imple¬ 
mented as functions acting on the event structure 
reading the incoming particles and producing the 
temporary ones. The type of the interaction is 
chosen according to the ratio of the total cross 
sections. 

In our presentation we focus on single pion 
channels. We present many comparisons with 
the existing experimental data. In the near fu¬ 
ture the generator will be supplemented with a 
module with nuclear effects. 

2. FRAGMENTATION ALGORITHM 

In our MC we use the DIS formalism to gener¬ 
ate events in the whole kinematical region where 
inelastic reactions are possible. In order to get 
the event record for the final state we assume 
that interaction occurs always on a particular 
parton and then the fragmentation of interact¬ 
ing quark and spectator is performed by means 
of PYTHIA6 |H1 routines. 
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Figure 1. Charged hadrons multiplicities for 
vp —> pL~. Data points taken from |H] are 
connected by lines. The results of our simulations 
are shown as separated points. 


Figure 2. Comparison of CC 1-pion functions 
with elasticities of N and A resonances. For neu¬ 
tron the sum of 1-pion functions for two exclusive 
channels is shown. 


The inclusive cross section for the scattering off 
nucleon is given by 
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Structure functions are assumed to be those de¬ 
fined in the parton model i.e. the combinations 
of a PDFs 

Fi {x, [qj {x, Q^) + q, (x, Q^)] 

3 

^3 {x, = 2 ^ [qj (x, - qj (x, Q^)] (3) 

3 

F 2 {x,Q^) = 2xFi{x,Q^) 

Using structure functions the cross section @ 
is rewritten in terms of contributions from sepa¬ 
rate partons qi 




(4) 


where Ki is a kinematic factor for parton qi. 


A probability of reaction on a given quark is: 


d3a’i'/dxdy 
- J2d^(7^^/dxdy 


( 5 ) 


PYTHIA fragmentation routines require a sys¬ 
tem of quark and diquark and perform fragmen¬ 
tation and hadronization using the LUND algo¬ 
rithm. Depending on the interacting parton, we 
distinguish several cases [ 7 ]: 


• In the case of the scattering off the valence 
quark, a string is formed from the created 
quark and the remaining diquark. 


• In the case of scattering off a sea quark u 
or d, the remaining anti-quark annihilates 
with appropriate valence quark, and the 
created quark forms a string with the the 
remaining diquark, exactly as in the previ¬ 
ous case. 


• If scattering off an anti-quark u gives an 
anti-quark d or scattering off an anti-quark 
d gives an anti-quark u, the created parton 
annihilates with a valence quark. 


dxdy 
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Figure 3. CC cross section for vN ^~X inter¬ 
action. Total cross section is split into contribu¬ 
tions from quasi-elastic, SPP [W < 2 GeV) and 
more inelastic processes. Data points are taken 
from ca, na, na, ca, m 
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Figure 4. CC cross section for DN inter¬ 

action. Total cross section is split into contribu¬ 
tions from quasi-elastic, SPP [W < 2 GeV) and 
more inelastic processes. Data points are taken 

from 113) II2I) CHI) CHI) E3) E3 



Figure 5. Cross section for vp fi~TT'^p. For 
data points and for simulations only events with 
hadronic mass TV < 2 GeV were included. Data 
points are taken from CHI) CS!) m. CSl) CHI) 

m 


• If scattering off an anti-quark u gives a 
strange anti-quark s or scattering off an 
anti-quark d gives an anti-quark c, it cre¬ 
ates with one of valence quarks a strange or 
a charm meson and the remaining quarks 
form a string for the fragmentation. 

• In the cases of scattering off a strange quark 
or anti-quark, the remaining strange con¬ 
stituent creates a strange meson with one 
of valence quarks and the remaining quarks 
form a string for the fragmentation. 

We fine tuned the PYTHIA6 generator param¬ 
eters of the fragmentation. In fig. (1) the com¬ 
parison of the charged particles multiplicities 

P{rich) = o-{nch)l ^ a-irich) ( 6 ) 

as obtained from our simulation with the data 
from the Fermilab bubble chamber is shown. 

3. 1-PION FUNCTIONS 

The only resonance we consider is the A and 
we have to estimate the single pion production 
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Figure 6. Cross section for vn ^~TT'^n. For 
data points and for simulations only events with 
hadronic mass W < 2 GeV were included. Data 
points are taken from ca, ca, ca, m 


Figure 8. Cross section for NC SPP channels. For 
data points and for simulations only events with 
hadronic mass W < 2 GeV were included. Data 
points are taken from m 
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Figure 7. Cross section for un —> /r“7r°p. For 
data points and for simulations only events with 
hadronic mass W < 2 GeV were included. Data 
points are taken from ca, d, m. m 


cross section as a fraction of the inclusive DIS 
cross section extrapolated into the resonance re¬ 
gion. This is done separately for each SPP chan¬ 
nel and the obtained fractions are called 1-pion 
functions. They are the probabilities that in a 
given point in the kinematically allowed region 
the final state is that of SPP. 

d^gPIS-SPP 

(7) 

dWdv 

In our generator are reconstructed using 

the LUND fragmentation algorithm. They turn 
out to be functions of W only and are shown in 
fig. (2). We see that up to the threshold for 
two pion production, 1-pion function for proton 
and the sum of functions for neutron are equal 1. 
In more common langauge 1-pion functions can 
be recognized as average elasticities of resonances 
^N-Tr/^total 0. In fig. (2) we see that in fact the 
values of 1-pion functions are close to resonance 
elasticities in a wide range of hadronic invariant 


mass. 
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Figure 9. Distribution of events in hadronic mass 
for BNL experiment [H] and predictions of our 
Monte Carlo for vp pL~'K^p 


Figure 10. Distribution of events in hadronic 
mass for BNL experiment and predictions of 
our Monte Carlo for vn —>■ p,~TT^n 


4 . SINGLE PION PRODUCTION 


Our model of SPP combines in a smooth way 
the A excitation model with the SPP part of 
the DIS cross section. We choose a linear tran¬ 
sition with respect to hadronic invariant mass 
W € (1.3,1.6) GeV. As a bonus we obtain an ar¬ 
tificial resonance-like behavior of the cross section 
at IF ^ 1.5 GeV which closely resembles the con¬ 
tribution from the Di 3 ,S'ii resonances ^D]. We 
describe the non-resonant background as a small 
admixture of the DIS SPP contributions at low 
values of W. Our MC reproduces the following 
analytical expression for the cross section: 


dW 


where 


da^ 

dW 


(1 


a(W)) 


da^^^ 

dW 


F^^^(W)a{W) 


( 8 ) 


W - Wth 

a{W) = 0(1.3GeF-VF)—- ^ao 

Irmm — yyth 

+ eiW^a. - IF)0(IF - IF™„) (9) 

W — Wmin + Oto{Wmax ~ W) 

+ 0(VF-IF,„a.) 


ao S (0, 0.3), depending on the channel. We note 
that a similar value for the division line in the 


hadronic invariant mass between resonance and 
DIS contributions = 1.5 ±0.02) 

was found by Naumov et al. cn by fitting proce¬ 
dure to the existing set of experimental data. 

The performance of our generator is presented 
on a series of plots. First we show contributions to 
the inclusive cross section for neutrino and anti¬ 
neutrino interaction on isoscalar target (figs. 3- 
4) from three theoretically separated dynamical 
mechanism. The SPP contribution is restricted 
by a cut IF < 2 GeV. 

The cross sections for GG SPP channels are 
shown in figs. 5-7. We conclude that the agree¬ 
ment with the data is satisfactory. In fig. 8 we 
show the plots for NG SPP channels. 

We also compared the distribution of events in 
hadronic mass for SPP channels with the data 
from the BNL experiment. It is an important test 
because our procedure of modelling SPP chan¬ 
nels is different from what is done in other MG 
codes. We used the BNL neutrino beam and gen¬ 
erated the same number of events as reported in 
CSl The results are shown in figs. (8-10). In 
the case of neutrino-proton reaction the agree¬ 
ment is excellent. In the case of vn p,~TT^p 
reaction the agreement is very good but our sim¬ 
ulations give too high A peak. In the case of 
vn —> /i“7r+n there is an experimentally mea¬ 
sured access of events with small invariant mass 
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Figure 11. Distribution of events in hadronic 
mass for BNL experiment m and predictions of 
our Monte Carlo for z/n ^ fj,~Tr^p 

(smaller then Ma = 1232 MeV) which is not re¬ 
produced by our simulations. 

5. FINAL REMARKS 

We find our results for SPP encouraging. An 
improvement in the vn —> p,~TY^n channel can 
probably be achieved by a more accurate treat¬ 
ment of the non-resonant background. 


Acknowledgments 

The collaboration of Krzysztof Graczyk at 

early stages of the MC generator and many useful 

discussions are acknowledged with pleasure. 

REFERENCES 

1. G. Battistoni et ah, http://nuint.ps.uci.edu/ 
proceedings/sala.pdf. 

2. C.H. Llewellyn Smith, Phys. Rep 3, no 5 
(1972) 261. 

3. E.A. Paschos, Ji-Young Yu, M. Sakuda, 
Phys.Rev. D69 (2004) 014013. 

4. M. Gluck et al. Z.Phys. C67 (1995) 433. 

5. A. Bodek, U.K. Yang, Nucl. Phys. (Proc. 
Suppl.) 112 (2002) 70. 

6. T. Sjostrand et al., Computer Phys. Com- 
mun. 135 (2001) 238. 


7. G. Ingelman, A. Edin, J. Rathsman, Gomput. 
Phys. Commun. 101 (1997) 108; F. Sartogo, 
PhD thesis (in Italian), supervisor P. Lipari, 
Rome 1994/95. 

8. D. Zieminska et al. Phys.Rev. D27 (1983) 47. 

9. O. Lalakulich, private communication. 

10. O. Lalakulich, http://users.ift.uni.wroc.pl/ 
~mb20 /presentations/Lalakulich.pdf 

11. K.S. Kuzmin, V.V. Lyubushkin, V.A. Nau¬ 
mov, How to sum contributions into the total 
charged-current neutrino nucleon cross sec¬ 
tion, hep-ph/0511308 

12. S.J. Barish et al. Phys. Lett. B66 (1977) 291; 
S.J. Barish et al. Phys.Rev. D19 (1979) 2521; 
G.M. Radecky et al. Phys.Rev. D25 (1982) 
1161; J. Gampbell et al. Phys.Rev.Lett. 30 
(1973)335, G.M. Radecky et al. Phys. Rev. 
D25 (1982) 1161 

13. P.G. Bosetti et al. Phys.Lett. B70 (1977) 273; 
D.C. Colley et al. Z. Phys. C2 (1979) 187; 
P. Bosetti et al. Phys.Lett. BllO (1982) 167; 

M. A. Parker et al. Nucl.Phys. B232 (1984) 1; 
P. Allen et al. Nucl.Phys. B264 (1986) 221; 
D. Allasia et al. Nucl.Phys. B343 (1990) 285 

14. C. Baltay et al. Phys.Rev.Lett. 44 (1980) 916; 

N. J. Baker et al. Phys.Rev.D25 (1982) 617 

15. T. Kitagaki et al. Phys. Rev.D34 (1986) 2554 

16. D. MacFarlane et al. Z. Phys. C26 (1984) 1; 
P.S. Auchincloss et al. Z.Phys. C48 (1990) 
411; W.G. Seligman et al. Nevis Report 
(1996) 292 

17. J.P. Berge et al. Z. Phys. C35 (1987) 443 

18. M. Jonker et al., Phys. Lett. B99 (1981) 265; 
J.V. Allaby et al. Z.Phys. C38 (1988) 403 

19. G.N. Taylor et al. Phys. Rev. Lett. 51 (1983) 
739; A.E. Asratian et al. Phys.Lett. B137 
(1984) 122; J. Bell et al., Phys.Rev.Lett. 41 
(1978)1008 

20. A.S. Vovenko et al., Sov. J. Nucl. Phys. 30 
(1979) 528. 

21. V.B. Anikeev et al., Z.Phys. CTO (1996) 39 

22. H.J. Grabosch et al., Z.Phys. C41 (1989) 527 

23. W. Krenz et al., Nucl.Phys. B135 (1978) 45 










